Large momentum part of fermions with large scattering length 
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It is well known that the momentum distribution of the two-component Fermi gas with large 
scattering length has a tail proportional to 1/fc 4 at large k. We show that the magnitude of this 
tail is equal to the adiabatic derivative of the energy with respect to the reciprocal of the scattering 
length, multiplied by a simple constant. This result holds at any temperature (as long as the 
effective interaction radius is negligible) and any large scattering length; it also applies to few-body 
cases. We then show some more connections between the 1/fc 4 tail and various physical quantities, 
in particular the rate of change of energy in a dynamic sweep of the inverse scattering length. 
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I. THE THEOREM 

Consider an arbitrary number of fermions with mass 
m in two spin states f and j, having an s-wave con- 
tact interaction with large scattering length a (|a| 3> ro 
0) between oppositite spin states, and confined by any 
external potential V ex t(r). To discuss their momentum 
distribution, it is convenient to use a box with volume f2 
and impose a periodic boundary condition. For a large 
uniform gas, fl is its actual volume; but for a gas con- 
fined in a trap, it is better to use (conceptually) a very 
large box which contains the region of the gas. Now con- 
sider any stationary state, that is, any energy level or any 
incoherent mixture of various energy levels (such as an 
equilibrium state at any temperature [2}). The momen- 
tum distribution (< 1 because of Pauli exclusion), 
that is, the average number of fermions with wave vector 
k and spin a, has in general a well-known tail of the form 
«ko- ~ C/k 4 at large k, where 

C = lim fc 4 nk(T 

k — >oo 

is independent of a [3]. In this equation, the limit is 
physically restricted to k <C 1/Vo, where ro is the effective 
interaction radius 

This tail determines how many large-momentum 
fermions (on average) we can find, if we measure the 
momentum distribution of the gas. This measurement 
is practically possible and has already been carried out 
to some extent 0, Q in some recent experiments on 
the BCS (Bardeen-Cooper-Schrieffer superfluid) to BEC 
(Bose-Einstein condensate) crossover behavior of such a 
Fermi gas |6(. The number of fermions with momenta 
larger than hK is, on average, 



N, 



(k>K) 



nc 



plus higher order corrections which are negligible at large 
K. Here K > l/l n , K > l/X dB , but K « l/r ; l n is 
the characteristic distance between two fermions in the 
opposite spin states 0, and \ c ib ~ ft/ \ffnkBTtemp is the 
thermal de Broglie wave length. H is Planck's constant 



divided by 2w, k B is Boltzmann's constant, and Ttomp is 
temperature Q]. 

If we tune the reciprocal of the scattering length @ 
adiabatically (very slowly), but keep the confinement po- 
tential [such as a trap or an optical lattice (for ultracold 
fermionic atoms)] fixed, the total energy E (the sum of 
the kinetic energy, the interaction energy, and the exter- 
nal trapping energy) of the gas will change accordingly, 
and we have a quantity dE/d(l/a), which we call the adi- 
abatic derivative of energy (with respect to the inverse 
scattering length 1/a). 

The above two quantities, the magnitude of the 1/k 4 
tail, and the adiabatic derivative of energy, can be inde- 
pendently measured in experiment. 

We show in this paper that these two quantities have 
a simple, exact and universal relation: 



dE 



h 2 nc 



Atttti d(-l/a)' 



(1) 



and we shall call it adiabatic sweep theorem. 

This theorem holds whenever the effective interaction 
radius ro l] is negligible, compared to the other rele- 
vant length scales of the problem: the scattering length, 
the average interparticle spacing, the thermal de Broglie 
wave length, and the characteristic length scale over 
which the external potential is inhomogeneous. In this 
limit (ro — > 0), the interaction between two fermions in 
opposite spin states is s-wave contact interaction, char- 
acterized by the scattering length only, and there is no 
interaction between fermions in the same spin state be- 
cause of Pauli blocking. 

The theorem applies to any stationary state, such as 
the ground state, any excited state, and any finite tem- 
perature state [J in thermal equilibrium. Generally 
speaking, if !</>,) are a set of energy levels of the system, 
then any mixed state described by a density operator 
p = pi\4>i) {4>i\ satisfies this theorem. The numbers of 
fermions in the two spin states iVj and the external 
potential (as long as it is fixed), and the scattering length 
are all arbitrary, provided that the conditions specified 
in the previous paragraph are satisfied. 



2 



II. THE PROOF 

Although the theorem can be proved with conventional 
means as well, we shall use the formalism developed in 
This formalism, as we shall see below, allows us to 
prove this universal theorem in a universal (and simple) 
way. 

The basis of the formalism developed in Q is a pair 
of linearly independent generalized functions A(k) and 
L(k), satisfying 

A(k) =1 (|k|<oc), =0, (2a) 

L(k) =0 (|k|<oo), f$F±jP =1, (2b) 



A(-k) =A(k), 



L(-k) =£(k).(2c) 



There is no contradiction at all, because the integrals of 
these functions over the whole momentum space are not 
solely determined by their values at finite k's. Like the 
delta function S(x), which has a singularity at x = 0, the 
above functions have appropriate singularities at infinite 
momenta such that the above integrals hold. 

This formalism is equivalent to the one based on the 
Fermi pseudopotential 4:nh 2 aS(r)(d/dr)r/m [?1 ITflL Hlj. 
but it has a greater simplicity and flexibility, making 
concrete calculations more straightforward, and allow- 
ing for some useful formal derivations. The present au- 
thor has used this new formalism to find another ex- 
act theorem and solve the two-body and three-body 
problems Q. The dimer-fermion scattering length cal- 
culated by the author 3] is completely consistent with 
independent calculations by others 0, 0], but has a 
much higher accuracy |3j. In addition, the author used 
it to calculate the dimer-dimer scattering length (unpub- 
lished), and got add = 0.6a, consistent with [ljj. The 
author also used this formalism to study the three-boson 
problem with large scattering length, the beyond-mean- 
field ground state of the dilute Bose gas, and the ground 
state of the dilute two-component Fermi gas with short 
range interaction. In all the instances, the results ob- 
tained within this seemingly paradoxical formalism are 
perfectly consistent with the celebrated results in the lit- 
erature. Finally, the formalism and the theorem derived 
in have been easily generalized to two dimensions . 

In the language of these generalized functions, the sys- 
tem's hamiltonian is 



H(a) = T + aU + V, 
where the three operators 
T = J2(^k 2 /2m)cl a c ka , 



(3) 



U = 



ka 



C q/2+kT C q/2-kl C q/2-k'i C q/2+k'T A ( k ')' 



qkk' 



V = Y d f d 3 r^ xt (r)V4(r)^(r) 



are all independent of a. Here c k(T and Vv( r ) are the con- 
ventional fermionic annihilation operators in momentum 
and coordinate representations, respectively. 

Any energy level |</>) associated with the hamiltonian 
H(a), such that H{d)\4>) = E\<fi), automatically satisfies 
a short-range boundary condition 3] 



q/2-kl C q/2+kTl^) — ^ ( an y 



(4) 



where ^(k) = A(k) + L(k) /Ana. This is just a rephrasing 
0, 0] of the familiar statement that the wave function is 
of the form coeff.(l/r — 1/a) + 0(r) when the distance r 
between two fermions in opposite spin states approaches 
zero. 

Further, it is shown that any state \4>) which satisfies 
Eq. must satisfy an energy theorem 0]: 



\H{a)\4>) = {(t>\{TiT + V)\4>), 



(5) 



for 



where wT is the shorthand 
Ek CT ??(kX^72m)cLv 

If we thoroughly examine the logic in the long proof 
of this energy theorem in [3j, we find that the symbol n 
is only associated with the boundary condition satisfied 
by (<f>\, but has nothing to do with \<p); that is, we can 
consider any energy eigenstate \<fi) of H(a) [such that 
H(a)\4>) = E\<j>)] and any energy eigenstate \<j>') of an- 
other hamiltonian H(a') [such that H{a')\<t>') = E'\cf>')}, 
where H[a!) differs from H(a) only in the scattering 
length, and prove that 



\H{a')\<t>') = {4>\{vT + V)W) 



(6) 



using the same logic as in the proof of the energy theorem. 

We may also follow another approach in 
0, and derive this equality much more easily: 
(</>[ (T + a'U + V)\4?) = {4>\(r)T + a'nU + V)\4?) = 
(<f>\(i]T + V)\4>'), where rjU is the shorthand for 

"raff Eqkk' ^( k ) C q/2+kT C q/2-kl C q/2-k'l C q/2+k'T A ( k/ )' 

and (4>\rjU = because of the boundary condition 

(01 Ek 7 7( k ) C q/2+k T C q/2-kl = - 

Now Eq. © can be further simplified as follows, be- 
cause \4>') is an eigenstate of H(a'): 



£W) = (0I(»7T + W>- 
Similarly, we have proved that 

£(0'|0) = (<f/\tfT + V)\<l>), 



(7) 



(8) 



where n' (k) = A(k) + L(k)/47ra'. 

Taking the complex conjugate of Eq. ©, noting that 
both rj'T and V are hermitian, and subtracting the result 
fromEq. £7J, we get {E'-E){<j>\tf) = ((f)\(r]-r]')T\^), and 
both the operator V and the function A(k) are canceled. 
What remains is 

{E'-E){ct>W) 
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h 2 k 2 



= (1/4™ - l/Awa'M ]T LM—^c^W). (9) 



Now we apply Eq. Q to a special situation, in which 
we do an adiabatic sweep of the inverse scattering length 
from 1/a to 1/a', and the state \<f>) is continuously trans- 
formed to (that is, |</>') is adiabatically connected to 
|0)). In the limit 1/a' — » 1/a, such that |0') — > |</>), we 
obtain a differential relation 



d£ 



d(l/a) 87rm 



$>(k)fc 2 



"k<r* 1 



By the basic property of L(k) 0, Ek £ ( k )/( k ) = 
fi/[d 3 A;/(27r) 3 ]L(k)/(k) = Q lim^^ /c 2 /(k) for any 
conventional function /(k). So the summation in 
Eq. ljTU|) is equal to Q um k— >oo 
Therefore, 



2f7C. 



d# 



n 2 nc 

47TTO ' 



and we have proved the adiabatic sweep theorem in the 
case of any pure energy eigenstate. 

In any incoherent mixed state of the form p = 
Si Pi 1 4>i)(4>i I ? where \<j>i) are energy eigenstates, both E 
and C are statistical averages of the E and C values in 
the states \4>i) with weights pi, so the adiabatic sweep 
theorem obviously still holds. 



III. TWO-BODY CASE: A CONVENTIONAL 
PROOF 

It is cumbersome to prove the above theorem with con- 
ventional means. Here we just demonstrate the two-body 
case, and then sketch the extension to the iV-body cases. 
This shall help to convince the readers that the adiabatic 
sweep theorem is indeed correct, even if they are not used 
to the succinct (and "unconventional" ) proof in the last 
section. 

Now consider two fermions in opposite spin states with 
an external potential T4 x t(r). We consider an energy 
eigenstate described by a wave function 0(ri,r 2 ), with 
scattering length a and energy E; we also consider an- 
other energy eigenstate described by a wave function 
$'(ri,r 2 ), with scattering length a' and energy E' . Here 
ri is the position of a spin up fcrmion, and r 2 , spin down. 



h 2 v\ h 2 v 2 2 . . , / 

2 +V ext (ri) + V ext {r 2 ) 



2m 



2 m 



4> = E4> (li) 



whenever ri ^ r 2 . There is also a short-range boundary 
condition: 

0(ri,r 2 )=A(r o )(l/r- l/o) + O(r) (r -» 0), (12) 



where ro = (ri +r 2 )/2 and r = ri — r 2 . Similar equations 
hold for <h' . It then follows that 



2m 



Witf*) = (E-E')^*<f>. (13) 



Making the substitution Vi = Vo/2 + V r and V2 = 
Vo/2 — V r (where Vo and V r are partial differential op- 
erators associated with r and r, respectively), we get 



c k J0)/<</#). (io) -JL- 



4m 
4V 



Vo-(0'*V o 0-0Vo0'*) 



0V r 0'*) =(E-E')cf>'*<j>. (14) 



We then integrate this equation over V t (e > 0), the set 
of all configurations (ri, r 2 ) in which r > e. The left side 
is converted to surface-flux according to Stokes theorem. 
Since the wave functions either decay to zero at infinity 
or satisfy a periodic boundary condition in a box, only 
the flux on the hypersurface r = e remains. Using the 
short-range boundary conditions, such flux is found to be 
(47rfi 2 /m)(l/a' — 1/a) J d 3 r A'*(r )^4(ro) plus correction 
terms which vanish in the limit e — > + . Taking the limit 
e — > + , we get 

[E - E') J d 3 rid 3 r 2 0'*0 

= (4^ 2 /m)(l/a' - 1/a) / d 3 r A'*(r )A(r ). (15) 



Now we specialize to the case in which <fi' is adiabatically 
connected to <fi, and take the limit 1/a' — * 1/a; we further 
assume that (f> is normalized. The above equation is then 
simplified as 



dE 



47rft 2 d(l/a) 



= J d 3 r |A(r )| ; 



(16) 



Now we turn to the evaluation of flC in the quantum 
state <f>. We need to transform the wave function to mo- 
mentum space: 



9 k!k 2 



d 3 rid 3 r 2 



(rir 2 )exp(-ikiTi-ik 2 T 2 ), (17) 



and the average number of spin-up fermions with wave 
vector k is 



k' 



>kk' 



(18) 



At large k, 0kk' is dominated by the contribution from 
the most singular term in 0(rir 2 ) (it is a general property 
of Fourier transformation), namely A(vo)/r. So 



W = fl^Aik + k')167r/(k - k' 



(19) 



plus higher order terms, at large k. Here A(q) = 
J d 3 ro^4(ro) exp(— iq • ro). Note that A(ro) describes the 
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center of mass motion of the two fermions (when their 
distance is zero) and should be everywhere finite and con- 
tinuous; consequently, A(q) should decay sufficiently fast 
at large q. It then follows that the summation in Eq. Ijl8|l 
should be dominated by terms in which |k + k'| <C k, at 
large k. So 



n kT = 53["~ 1 ^(q)V* a ] : 



(20) 



plus higher order corrections, at large k. So 

C= lim fc 4 n kT = fr^TT 2 / d 3 r |^(r )| 2 , (21) 

k^oo J 

and the same formula results if we study the spin-down 
fermion. Comparing this with Eq. (|16|l . we find 



h 2 nc 



AE 



If the two fermions are in the same spin state, there is 
no interaction between them, and E is independent of a; 
meanwhile, C — 0. 

The conventional proof of Eq. Q for more fermions is 
similar to that of the two-body case. The major addi- 
tional considerations are: 

• In the intermediate steps, the set T> € (see above) 
also excludes the configurations in which three or 
more fermions cluster in a tiny region of size e. 
When we apply Stokes theorem, we also need to 
consider surface contributions from these excluded 
configurations. But Pauli blocking suppresses the 
probability amplitude for three or more fermions 
to cluster in such a tiny region. It then follows 
that surface contributions from these three-body 
and more-body clusters vanish in the limit e — * 0, 
leaving us with contributions from the two-body 
clusters only. An equation similar to Eq. I|15|) thus 
results. 

• Three-body or more-body clusters can only con- 
tribute a tail of the form 0(1/ k a ) (a > 5) to the 
momentum distribution at large k, because of Pauli 
blocking 16] . The 1/fc 4 tail is solely due to the 
two-body (opposite spin) clusters, and an equation 
similar to Eq. l(2T|l still holds. 



IV. THE IMPLICATIONS 

A. Physical meaning 

The meaning of Eq. Q is intuitively clear. If we tune 
the inverse scattering length slightly and slowly, quan- 
tum mechanical first-order perturbation determines the 
amount by which the energy level shifts. In the limit of 
zero effective interaction radius, two fermions (in oppo- 
site spin states) do not interact unless they appear at 



the same position, so the energy shift should be propor- 
tional to the probability that this occurs. Meanwhile, 
f2C characterizes this probability recall that when 
two fermions appear at the same place, their momenta 
are large and nearly opposite. The definitive form of 
Eq. is of course not obvious. 



B. Large momentum part in the BEC-BCS 
crossover 

In this subsection we are restricted to a special (but 
interesting) case: the uniform Fermi gas with equal pop- 
ulations of the two spin states at zero temperature. The 
scattering length a is large (|a| 3> r ) and arbitrary. 

There is only one independent dimcnsionless parame- 
ter, x = — l/k F a, where 



kp — (37r 2 n 



iV3 



is the Fermi wave vector, and n is the number density 
of fermions. When x is negative and large, the ground 
state 17] is a BEC of fermionic dimers ]g; when x is 
positive and large, the gas becomes a BCS superfiuid |(|; 
in between, the system is in a strongly correlated many- 
body state, but there is no quantum phase transition . 
For all finite x values (including x = 0) , the ground state 
is believed to be a superfiuid 

According to Eq. QJ, we just need to take the deriva- 
tive of the ground state energy with respect to —1/a, 
to get the magnitude of the large momentum part, the 
contact intensity C. 

BEC limit- the ground state energy is [HI El GJ, El 



E/Vi = n d E d + 2-Kh 2 n d a d dlm d 

x (1 + 4.8144^/5 + 19.6539gln g + c 3 g 



(22) 



where n d = n/2 is the dimers' number density, E d — 
—H 2 /ma 2 is the energy of an isolated dimer at rest, 
m d = 2m is the dimer's mass, g — n d a dd is the Bose 
gas parameter. add = c P a Il4l is the dimer-dimer scat- 
tering length, and c P w 0.6 [lj|. The constant C3 is still 
unknown, since it depends on the difficult three-dimer 
scattering problem [20I . The contact intensity thus 
has an expansion 



C = Ann/a + c P ir 2 n 2 a 2 [l + 12.0360^ 

+ 78.6156gln.g + (4c 3 + 58.9617)g - 



(23) 



Unitary regime l/k F a is small or zero (the latter 

case is called unitarity limit 1231 l24l |25| ) . There is an 
expansion [H H |H El El llj 



E 

n 



5 2m 



k F a 



5v 



3(k F a) 2 



(24) 



where the current theoretical estimate of the constants 
are £ w 0.44 [H, |27|, C « 1 El IH and v « 1 
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26, 28]. Consequently, 



C = kU 2C/5tt 



4u/3tt 
k F a 



(25) 



(26) 



(27) 



is a universal constant. If ( w 1 |2fl |2JJ |28j , j « 0.13; 
this value of j is consistent with a recent Monte Carlo 
calculation of the full momentum distribution |2|j. Of 
course, if £ is appreciably different from 1, j must be 
appreciably different from 0.13. 

BCS limit- the ground state energy is [23| 



In particular, in the unitarity limit, 

C — jk F (unitarity limit), 

where 

J = 2C/57T 



E 



_3_ 

10 



1 



3?r 

o.oo9i4(fc F a) 3 



fcpa + 0.055661 (fc F a) 2 



0.018604(fc F a) 4 + 



(28) 



where 2(11 - 21n2)/357r 2 w 0.055661 is exact 111, but 
the last two coefficients still have uncertainties |22|; the 
BCS superfluid energy saving is exponentially small [~ 
exp(— 7r/2fc F |a|)] (see for example |30|'l and does not show 
up in this asymptotic power series of k F a |2^|. So 



C = 4ir 2 n 2 a 2 



1 + 1.049190fc F a 



0.2584(fc F a) 2 - 0.70135(/c F a) 3 



(29) 



where 12(11 - 21n2)/357r » 1.049190. 

The full momentum distribution in the BCS limit has 
actually been calculated by Belyakov [3l| and later cor- 
rected by Sartor et al [3^ | , to the lowest nontrivial order 
[0(a?)\ in the scattering length only, although these au- 
thors studied the problem in the context of a dilute hard 
sphere Fermi gas (with positive scattering length), and 
they did not foresee the applicability of their result to a 
dilute (l/k F ^> \a\) Fermi gas with a large (\a\ ^> r ) and 
negative scattering length |33| . The universal momentum 
distribution function [to 0{a 2 )\ they [3lll3^ | derived has 
a tail consistent with the leading order term in Eq. (|29|l . 

Note that C oc n 2 a 2 in the BCS limit, in contrast with 
the predictions based on the BCS equations, which lead 
to an exponentially small value [~ exp(— 7r/fcp|a|)] of C. 
BCS equations do not predict the correct behavior of C 
in the BCS limit. This observation is not incompatible 
with a recent experiment [5j , in which the measured large 
momentum tail of a trapped Fermi gas in the BCS regime 
is found 5] to be stronger than the theoretical prediction 
based on the BCS theory. 

These various expansions of C are plotted in Fig.^ It 
is almost certain that C monotonically decreases from the 




FIG. 1: (color online) Intensity of the large momentum part 
C = linik_ 00 k 4 n ktJ in the BEC-BCS crossover at zero tem- 
perature. BEC reg ime: Eq. I|23|l . in which the dashed line 
is for c 3 = 140 |20L solid line for c 3 = 40 and gray 
dotted line, c 3 = 30 |22| (the upward bending of each of these 
curves toward the unitary regime marks the breakdown of the 
low-density expansion). Unitary regime: Eq. I25H combined 
with the numbers £ fa 1 and Dfsl (note that the actual slope 
of the curve in the unitarity limit is probably less steep; see 
text). BCS regime: Eq. 129L which should be less reliable 
when —l/k F a w 1. 



BEC limit to the BCS limit, at a fixed particle density. 
The parameter v determines the slope of the C curve in 
the unitarity limit (—l/k F a = 0), and it seems that the 
value v « 1 (responsible for the slope of the short line in 
the unitary regime) is probably too large. 

What we have learned here, which is based on Eq. QJ, 
is remarkable: simply measuring the coefficient of the 
l/k A tail at various scattering lengths, one could pin 
down the whole equation of state of this novel Fermi gas 
0. 

Even the famous constant [see Eq. (|24|l ] 

characteri zing the ground state energy of the unitary 
Fermi gas |23t |24| , can be alternatively determined from 
measurements of C away from the unitarity limit. Sup- 
pose that 

C = k A F c{-l/k F a), 
and we can show from Eq. Q that 



c(x)dx. 



(30) 



c(x) also must satisfy a constraint [according to Eq. 
and the behavior of the energy density in the BEC and 
BCS limits]: 



\c(x) + 4a;/37r]da; + 



+ 00 



c(a;)dx = — . (31) 

57T 
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We finally remark on the practical accessibility of the 
C /fc 4 tail in the momentum distribution. In the BCS 
limit, the result in j^J indicates that the difference 
between k n^a and C is already as small as a few per- 
cent, when k/k F = 5, and it decreases fast at larger k 
values. The situation is similar in the unitary reg ime, as 
indicated by a recent Monte-Carlo simulation [2!j. In the 
BEC limit, we expect the difference between fc 4 nko- and 
the true value of C to be negligible when k is more than a 
few times 1/a (a characterizes the size of each fermionic 
dimer). 



C. Large momentum part of the trapped gas 

For the gas confined in a trap, the parameter C in Sec.^ 
is only the average contact intensity over the arbitrary 
volume fi. It is the integrated contact intensity, S7C, that 
is uniquely defined. 

We can, however, uniquely define the local contact in- 
tensity 

C(r) = -4tt lim S V s 2 (V4(r)^(r + s)), (32) 

which is also independent of a. The limit is restricted to 
s»fO' It is straightforward to show that 

nC = J d 3 rC(r). (33) 

An advantage of C(r) is that for a trapped cloud of 
many fermions, we can use the local density approxima- 
tion to predict the value of CIC from the equation of state 
of the uniform gas. Let us illustrate this with the zero 
temperature case. First, the local density of the gas is 
estimated from the equation /iiocai[a(r)]+V^ x t(r) =const, 
where ii\ oca \(n) is the chemical potential of the uniform 
gas with density n. Then C(r) is approximately calcu- 
lated from this local value of n using the formula of C 
for the uniform gas. Finally we do the spatial integral 
and predict (approximately) the intensity of the large 
momentum part of the trapped gas. Note that in this 
analysis, the formula for C of the uniform gas is also 
solely determined by the equation of state. 

By comparing our prediction with the measured large 
momentum part, we could find out whether a particular 
equation of state we have used in this kind of analysis is 
correct and, if not, we could adjust it. In this way, the 
experimental data on the large momentum part may also 
help us to pin down the equation of state of the uniform 
gas. 

We now turn to another question: how to directly test 
the adiabatic sweep theorem in a trap at any tempera- 
ture 0? At least one approach is feasible in principle (for 
1/a < only) and involves almost no theoretical approx- 
imations. 1) We first prepare a sample of equilibrated 
Fermi gas in a particular trap and then take in situ im- 
age of the gas to determine its external trapping energy. 



2) We produce the trapped sample again, following ex- 
actly the same recipe as before, but then suddenly switch 
off the trap (without changing the scattering length) to 
allow the gas to expand; after enough time, an image 
is taken, and the kinetic plus interaction energy of the 
gas is deduced from the second moment of the spatial 
profile. The sum of the above two energies is the to- 
tal energy. 3) The sample is prepared the third time, 
and then the inverse scattering length is adiabatically 
tuned (over a time scale longer than the relaxation time) 
to a somewhat different value 1/a', and then the exter- 
nal trapping energy is measured in the same way as step 
one. 4) The sample is prepared and the inverse scattering 
length adiabatically tuned to 1/a' again, and the kinetic 
plus interaction energy is measured like in step two. The 
sum of the energies from steps three and four is the en- 
ergy of the sample at that new scattering length, which 
should be different from the sum of energies from steps 
one and two. This difference, according to our theorem, 
must determine the intensity of the 1/fc 4 tail (if the new 
inverse scattering length does not differ too much from 
the old one). So we can prepare the sample for the fifth 
time, measure its momentum distribution, and compare 
the large momentum part with the prediction based on 
our theorem. 

The above five-step approach (if somewhat tedius) 
should be robust. A variant of this approach (which may 
be more practical) is to measure the total energy at a 
series of —1/a values (all of the states are adiabatically 
connected), do a curve fit, and take the derivative of the 
curve with respect to —1/a. In this way, the adiabatic 
derivative of energy can be measured more accurately. 
No local density approximation is involved in this ap- 
proach. 

D. Thermodynamics at any temperature 

In this subsection we consider the uniform gas at any 
temperature |2|. The two spin states are assumed to 
be equally populated. Now we have two independent 
dimensionless variables, x = —l/k F a, and the average 
entropy per fermion s = S/k B N. Here S is the entropy, 
and N is the number of fermions. 

The contact intensity is then determined by a universal 
dimensionless function c(x, s): 

C = k 4 F c(x,s). (34) 

This function contains all the information about the ther- 
modynamic properties of this novel Fermi gas, including 
the phase transition between the superfluid and normal 
states. Here we use entropy per particle instead of the 
reduced temperature (temperature divided by the Fermi 
temperature), because entropy is invariant in an adia- 
batic sweep and is thus a more convenient variable here. 

The properties of this Fermi gas are known in each 
of these three limits: x — » — oo, x — > +oo, and s — > 
oo [35| . The only difficult region is at x ~ 0(1) and 
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s ~ O(l), abbreviated below as D. If we start from the 
x — > ±00 regions, and tune the inverse scattering length 
adiabatically into D, we know the entropy of the final 
state from that of the well-known initial state. Measuring 
C at various x values in between, we can integrate C over 
x to get the energy of the final state. If we repeat the 
same measurements at various initial entropies, we can 
calculate the derivative of energy with respect to entropy 
at any scattering length, and obtain the temperature of 
the gas in D. 

This program is viable in principle, and it has the 
advantage of being model-independent (it is needless to 
say that there are other experimental approaches to the 
thermodynamic quantities as well; purely for the sake of 
space, we will not list or compare them in this paper), 
since almost no theoretical approximations are involved 
j36| . In reality it is complicated by the fact that samples 
of such Fermi gas studied to date are all in small traps 
or in optical lattices, which have inhomogeneous density 
distribution. However using some ingenious data analysis 
method, we might still be able to extract valuable infor- 
mation from experimental data on the large momentum 
part. 

We can also derive a useful relation about the pressure 
of the uniform Fermi gas with large scattering length and 
equal populations of the two spin states, in thermal equi- 
librium: 



p= 2 h 2 c 

3^ E Ylitarn 



(35) 



where p E = E/VL is the energy density. This relation 
holds at any temperature Q. 

Proof. On dimensional grounds, the energy density is 
p E = H k F f(x,s)/m, where x — —l/k F a and s is the 
average entropy per fermion. Using Eq. lJT]l. we can eas- 
ily show that c(x,s) = 4n(d/dx)f{x,s), where c(x,s) is 
defined by Eq. pt|). 

Suppose that the total number of fermions is N, and 
the volume is fl. Now we compress the gas adiabatically 
(but keep the scattering length unchanged) , to reduce its 
volume to a slightly smaller value £7' = (1 — o) 3 SI. Con- 
sequently, the average interfermionic spacing l/k F is re- 
duced by a factor 1 — o, so x' = (1 — o)x. The Fermi wave 
vector becomes k F = k F /(l — o). So the energy of the gas 
changes by AE = n'h 2 kff(x', s)/m - m 2 k 5 F f(x, s)/m, 
and the pressure is P = — lim _>o AE/(Q! — fi) = 
(H 2 k F /m)[2f/3 - x(d/dx)f/3] = (h 2 k F /m)(2f/3 - 
xcjYlit) = 2p E /3 + h 2 C/(l2nam). □ 

If we further express the energy density in terms of the 
momentum distribution namely p E = h 2 C / Airam + 
J2 a f[d 3 k/(2TT) 3 }A(k)(h 2 k 2 /2m)n ka , we see that the 
pressure (in the equilibrium state) is a simple functional 
of the momentum distribution 



P 



h 2 C 
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A-Kam ^ K ^°° J\k\<K i^) 2 2m 



d 3 k h 2 k 2 
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h 2 c 
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d 3 k kl H 2 k 2 

A(kj— — n ko 



E 



d 3 /t 



;A(k) 



(27r) 2 "~ v "' 2m 
L(k)l h 2 k 2 



Aira 



2m 



(36) 



where K <C 1/V physically. We shall call this result 
pressure theorem (of the two-component Fermi gas with 
large scattering length). 



E. Dynamic sweep 

What if we tune the inverse scattering length at a fast 
rate, so that the system's evolution is not adiabatic? Do 
we still have some insight into the rate of change of the 
system's energy (or, more precisely, expectation value of 
energy, since the system is not in any energy eigenstate 
if it is changing with time) [37j? 

If the characteristic time scale i S woep of the sweep is 
comparable to or shorter than mr^jh |38| . we have to 
either modify or completely give up the s-wave contact 
interaction model; this topic will not be discussed here. 

If ^sweep mr^/H [33, one can reasonably expect 
that the sweep can still be regarded as adiabatic from 
the perspective of the short-range boundary condition, 
ze, at any given time, the wave function is of the form 
cocff.[l/r — l/a(t)] + 0(r) when the distance r between 
two fermions in opposite spin states is small. Here a(t) 
is the instantaneous scattering length. But since mr^jh 
is much shorter than the other characteristic time scales 



of the system [such as the Fermi time m/(hk F )], t sweep 
may still be comparable to or even shorter than those 
other time scales, so that the whole quantum state may 
not evolve adiabatically. 

Given that the above instantaneous short-range 
boundary condition is satisfied, we can show that the 
energy expectation value of the system changes at a rate 



dE _ h 2 nc(t) di-ait)- 1 } 

~dt ~ 



dt 



Jd 3 r Pa (rt) 



dKxt(rt) 

dt ' 

(37) 

where C(t) is the instantaneous coefficient of the 1/fc 4 
tail of the momentum distribution, p a (Yt) is the spatial 
density of fermions in the spin state a, and for general- 
ity we also allow the external potential to change with 
time. We assume that the inverse scattering length is a 
real number at any instant of time. We shall call Eq. I|37|) 
dynamic sweep theorem. This theorem holds for any pop- 
ulations (no matter equal or unequal, small or large) of 
the two spin states, any nonzero scattering length, any 
quantum state, and any sweep rate, provided that the 
s-wave zero-range interaction model is still valid [3|j . 

Proof. Let us first describe the time evolution of the 
quantum state. In the Schrodinger picture, a quantum 
state \4>(t)) satisfies the Schrodinger equation 

d\<f>(t))/dt = -(i/h)H(t)\ct>(t)) 
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and the instantaneous short-range boundary condition 



C q/2-ki C q/2+kT 



\4>{t)) = (any q), 



where H(t) = T+a{t)U+V(t) [T, U and V(t) are defined 
in Sec. m except that now V ex t can change with time], 
and r)(kt) = A(k) + L(k)/47ra(t). 

We then prove the equation for the conservation of 
probability: 



dt 



o. 



d<0(t)|0(*)>/dt = <0(t)|d|0(*)>/d* + c.c. - 

(-i/n)<0(t)|ir(t)|^(i)) + c.c, but (4>{t)\H{t)\<f>(t)) = 

(4>(t)\[r](t)T+V(i)]\^)(t)) according to the energy theorem 
yj , and it must be real since the instantaneous scattering 
length is assumed to be real. So d(</>(i)|</>(i))/di = 0. 

Given the probability conservation, we may choose the 
normalized quantum state, such that (<f>(t)\(f>(t)) = 1. 
Then the system's energy expectation value is simply 

25(f) = m)\H(t)\m- 

Now consider two instants of time, t and t' . For con- 
ciseness, we shall suppress the symbol t for a quantity 
at time t, and use a superscript ' for a quantity at t for 
distinction. We then have 



E' -E 



b'\H'\ 
V\H'\ 
MH'k 



,\H\4 
cj>\H\ 
b'\H\ 



R, 



where R = 



)Hl 



(4>\h'WY = mT + o!u + v')W)* 

= {<t>\(tfr + a'rjU + V')\4f)* 
= (4>\( V T + V')\ct>'r 

= {4>'\W + v')\<j>), 

and similary (4>'\H\<f>) = (<f>'\(r)'T + V)\<f>). So 

E'-E= (tf\[(n - rt)T + (V - V)] |0) + R, 

Finally we divide both sides of the above equation by 
t' — t, and take the limit t' — > t. We get 



dE 
~dt 



dr]T dV 
~dT + ~dJ 



dt 



but {d{4>\/dt)H\4>) = {i/K){H\<t>)yH\<t>) = (>lh)\H\4>)\ 2 is 
purely imaginary (or zero), so the two terms in the last 
bracket cancel. Equation l|37|l then easily follows (note 
that the fermionic annihilation and creation operators 
are independent of time in the Schrodinger picture). 

If the system is in a mixed state, it can be described 
by a density operator J2i Pi\4>i{t)) ( < Pi(t)\i where the pure 



states \4>i(t)) obey the Schrodinger equation and the in- 
stantaneous short-range boundary condition, YliPi = 1> 
and (4>i(t)\4>j(t)} — Sij (if this orthonormality condition 
is satisfied at a time, it will be satisfied at all other time, 
because of the probability conservation) . The three terms 
in Eq. H37fl are all statistical averages of the correspond- 
ing values in the pures states \<pi{t)) (with weights pi), so 
Eq. still holds. □ 

Note that in this proof, we never write expressions 
like H(ti)\4>(t2)) (t\ t 2 ), since such entities contain 
delta function singularity, and the inner products be- 
tween them and any usual state vectors are divergent. 
It is only when the scattering length in H coincides with 
the one in \<p) that the delta function singularity in H\<f>) 
vanishes, as is the case when H and \(f>) are for the same 
instant of time. 

Equation l|37() indicates that if V ex t is held invariant, 
but —l/a(t) decreases, the expectation value of energy 
also decreases with time. This appears to be incompat- 
ible with the common wisdom that if the sweep time is 
comparable to or much shorter than ma 2 /H, and a is pos- 
itive and decreases with time, a fermionic dimer can be 
dissociated and the energy can increase. The resolution 
to this paradox lies in the fact that in such a sweep, the 
dimer also has a chance to evolve into a more tightly 
bound molecule, and the expectation value of the energy 
decreases. For more details, see Appendix 1X1 

Inspired by Appendix^ we arrive at a general conjec- 
ture: if the scattering length is suddenly changed from 
a large value ao to another large value a\, over a time 
scale At (At is much smaller than ma\/h and the other 
time scales of the system except mr^/H), then OC(i) is 
nearly a constant during this interval, and approximately 
equal to its value right before the scattering length be- 
gins to change. This is valid for any number of fermions. 
We may exploit this feature, to measure the integrated 
contact intensity as follows. 

Consider any number of fermions in a trap having large 
scattering length ao and negligible r$. Suppose that the 
external trapping energy is E cxt , and the internal energy 
(kinetic plus interaction) is i?i n t . If we suddenly change 
the scattering length to a large and negative (or infinite) 
value oi, and rapidly turn off the trapping potential at 
the same time, the gas will expand, and we can determine 
its final energy E' from the second moment of the spatial 
profile after a long enough time. The dynamic sweep 
theorem, combined with the ansatz that QC(t) is nearly 
a constant in the short interval in which the scattering 
length changes, implies that 



E — Ew 



h 2 nc f 1 

47T771 V 



1 



a ax 



(38) 



where VLC is the integrated contact intensity right before 
the rapid change of the scattering length. So if we plot 
E' against l/a Q — I /ax, we will get a straight line, whose 
slope determines f2C, and whose vertical intercept equals 
the initial internal energy of the gas. It does not matter 
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whether cto is positive, negative, or infinite. But a\ had 
better be negative or infinite, or some dimers will survive 
the expansion of the gas, complicating the measurement 
of energy. 

Unlike another approach (depicted in Sec. HJ) to mea- 
sure VIC, this one makes use of the whole spatial profile 
of the gas after its expansion, and thus has a much lower 
noise-to-signal ratio. 

If At is very small, but not extremely so, it is possible 
to derive the lowest order correction to Eq. JSSJ. Here 
we just show the result; the complete proof will appear 
later. Suppose that 

b(t) = l/a(t) 

approaches a constant 60 sufficiently fast at t — * — oo, and 
approaches another constant b\ sufficiently fast at t — > 
+00; the time evolution of b(t) occurs around t = 0, over 
a characteristic time scale At, which is small compared 
to the other time scales of the system except mr\ j h. If 
rar\jh is negligible, 



AE = 



h 2 nc 



8h 



-Ab-\l— Vt 
nm J J 

t'<t 



V b(t')b(t)dt'dt 



(39) 

where AE = E' — Ei nt , Ab = 61 — 60, 6 is the derivative 
of 6 with respect to time, and Co is the value that C 
would have at t = if b had remained equal to bo for 
all t. Note that this lowest order correction to Eq. i|38|) 
scales like V~At, and is independent of the state of the 
fermions; the next order correction will most likely scale 
like At. In the derivation of this formula, one exploits 
the knowledge that C{t) only changes a little during the 
fast ramp of the inverse scattering length. Note also that 
this correction term to AE is not always negative; if, for 
example, 6 suddenly increases from 60 to 62, and after a 
short period At it suddenly decreases from 62 to 61 , then 
the correction term is positive. 

The author thanks T. L. Ho for a stimulating sugges- 
tion. The author thanks K. Levin for the introduction 
to ultracold Fermi gases, and S. Giorgini for communica- 
tion. 



APPENDIX A: DISSOCIATION OF A 
FERMIONIC DIMER AFTER A SUDDEN 
CHANGE OF THE SCATTERING LENGTH 

Our purpose here is to illustrate the dynamic sweep 
theorem with a simple two-body problem. 

A dimcr with fcrmionic scattering length oq > (and 
so » rj) is initially at rest. The initial value of OC is eas- 
ily calculated to be 8n/a n . At time t « 0, the scattering 
length is suddenly changed to a (\a\ 3> ro), over a time 
scale that is both much shorter than mmimag/h, ma 2 /H) 
and much longer than mr 2 , / h. What is the energy expec- 
tation value of the system afterwards? Suppose U cx t = 0. 



Answer: since C(t) can only change continuously with 
time (we will verify this below), it can be regarded as 
almost a constant during the short interval during which 
the inverse scattering length is changing. Equation 13 7|) 
can then be integrated to yield AE = h 2 HC(Q)(l/ao — 
l/o)/(47rm) = [h 2 /m)26 (6 - 6), where 

6 = l/a , b=l/a. 

The system's initial energy is —(tf/rnjb 2 ,. So its en- 
ergy expectation value after the change of the scattering 
length should be 



E = {h 2 /m){bl - 26 6) (OO). 



(Al) 



Now we confirm this result with the concrete calcula- 
tion. The system's wave function is 



^(r,t) 



= <p(r, t) 
\f2na~Q r 



(A2) 



which is normalized: f d 3 r|-0(r, t)| 2 = 1, and f(r,t) — 
exp(— 6 r + iHbot/m) for t < 0. Here r is the distance 
between the two constituent fermions. 

For t > 0, we can expand the wave function in terms 
of the new energy eigenstates (associated with the scat- 
tering length a), and get 

ip(r, t) — Aexp(— br + ihb 2 t/m) 

poo 

+ / B(k) sin[fcr — a(k)] exp (-ihk 2 t/m)dk, (A3) 
^0 

where the scattering phase shift a(k) = arccos , b ^ k 

6 ^ w ^/i+b 2 /k 2 

is in the range (0, 7r); this ensures that the phase shift 
does not jump when we pass through the unitarity 
limit. The base functions satisfy orthogonality condi- 
tions f exp(— br) sin[fcr — a(k)]dr = (if a > 0) and 
/ °° sin[fcr — a(k)] sin[fc'r — a(k')]dr — [tt /2)8{k — k') (any 
a), and k, k' > 0. 

The coefficients A and B(k) can be determined from 
the projection of tp(r, 0~) onto the new energy eigen- 
states, because the wave function can only evolve con- 
tinuously, ip(r, CU) = ip(r, + ). The result is 



.4 



B(k) 



260(6) 



•60' 
2(6 - b )k 



7r(fc 2 + b 2 )Vk 2 + b 2 '' 



(A4) 
(A5) 



where 9(b) = 1 if 6 > and 0(6) = if 6 < 0. 

These functions satisfy the probability conservation 
condition Pq + Pi = 1, where 

Po = A 2 6o/6 = 466 o 0(6)/(6 + 6 o ) 2 

is the probability for the dimer to remain in the bound 
state, and P x = J dP 1 = J °° 7rb [B(k)] 2 dk = (b 2 + b 2 - 
2|6|6o)/(6 + 6q) 2 is the probability for it to be dissociated. 



10 



The system's energy expectation value for t > is 
clearly 

E = (h 2 / m)^y—b 2 P a + J k 2 ■ Trb [B (k)} 2 dk^ . (A6) 

Substituting Pq and B{k) into this expression, we find 
E = (ti 2 /m)bo(bo — 2b), consistent with our previous ex- 



pectation based on the dynamic sweep theorem. 

Now we verify the continuity of the contact intensity. 

For any given t < 0, cp(r, t) — exp(ifi&ot/m)(l — b$r) + 
0(r) when r — ► 0. 

For any given t > 0, according to Eqs. (|A3|) . I|A4(1 and 

p(r,f) = v(0,i)(l - br) + O(r) when r -> 0, where 



¥>(0,*) 



6q exp(i6Qfe/m) < 1 — erf exp (m /A)boyJht/m > + b exp(ib 2 ht / m) < 1 + erf exp(i7r/4)6^/ ht/m 



ba + b 



r 



(t > 0). 
(AT) 



erf(,z) = (2/t/tt) J* exp(— w 2 )dw is the error function. 
The value of (p(0,t) (t > 0) in the special case b — 
—bo is given by the limit b — > — 6 . In any case, 
lim t _>o- </?(0,i) = lim t ^ + i/j(0, i) = 1, and this means 
that fiC(i) = (87r/ao)|y(0, i)| 2 is continuous. But note 
that [(d/dr)ip(r,t)] r= o is discontinuous at t = 0. 

The Fourier transform of <p(0,t) (— oo < i < oo) is 
simple: 



£(0,cj) = 



Po-P 



(u - ie + pg) (Vw + ie - ip) ' 
where p = yfi/m/ao, p = y/H/m/a ^ p , and e — > 



(A8) 



This two-body problem also helps to assure us that 
the momentum distribution of fermions with large scat- 
tering length (|a | S> r ) is indeed physically measurable 
|40j |. According to Eq. (|A5J) . i/ </ie dimer is dissociated, 
the momentum distribution of the fermions in the limit 
t — > +oo (corresponding to a directly-measurable spatial 
probability distribution) is 



1 



n kcr cx 



(l + a 2 fc 2 )(fc 2 + l/a 2 ) 



2\2 ' 



(A9) 



which in the limit \a\ <C ao approaches the momentum 
distribution of the fermions in the original dimer. 
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